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Possible Evidence For Dark Matter Annihilation In The Inner Milky Way From The 

Fermi Gamma Ray Space Telescope 
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We study the gamma rays observed by the Fermi Gamma Ray Space Telescope from the direc- 
tion of the Galactic Center and find that their angular distribution and energy spectrum are well 
described by a dark matter annihilation scenario. In particular, we find a good fit to the data for 
dark matter particles with a 25-30 GeV mass, an annihilation cross section of ~ 9 x 10~ 26 cm 3 /s, 
and that are distributed with a cusped halo profile, pir) oc r~ 1A , within the inner kiloparsec of the 
Galaxy. We cannot, however, exclude the possibility that these photons originate from an astro- 
physical source or sources with a similar morphology and spectral shape to those predicted in an 
annihilating dark matter scenario. 



Searches for dark matter annihilation products are 
among the most exciting missions of the Fermi Gamma 
Ray Space Telescope (FGST). In particular, the FGST 
collaboration hopes to observe and identify gamma 
rays from dark matter annihilations occuring cosmolog- 
ically [3, as well as within the Galactic Halo [2|, dwarf 
galaxies Q , microhalos 0] , and the inner region of the 
Milky Way @. 

Due to the very high densities of dark matter predicted 
to be present in the central region of our galaxy, the 
inner Milky Way is expected to be the single brightest 
source of dark matter annihilation radiation in the sky. 
This region is astrophysically rich and complex, however, 
making the task of separating dark matter annihilation 
products from backgrounds potentially challenging. In 
particular, the Galactic Center contains a 2.6 x 10 6 Mq 
black hole coincident with the radio source Sgr A* Q, 
the supernova remnant Sgr A East, and a wide variety 
of other notable astrophysical objects, including massive 
O and B type stars, and massive compact star clusters 
(Arches and Quintuplet). 

Since its launch in June of 2008, the Large Area Tele- 
scope (LAT) onboard the FGST has identified as pho- 
tons several hundred thousand events from within a few 
degrees around the Galactic Center. In addition to pos- 
sessing the effective area required to accumulate this very 
large number of events, the angular resolution and energy 
resolution of the FGST's LAT are considerably improved 
relative to those of its predecessor EGRET. As a result, 
this new data provides an opportunity to perform a pow- 
erful search for evidence of dark matter annihilation 0]. 

Dark matter annihilations are predicted to produce a 
distribution of gamma rays described by: 
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where < av > is the dark matter particle's self- 
annihilation cross section (multiplied by velocity), m^ m 
is the dark matter particle's mass, ip is the angle away 



from the direction of the Galactic Center that is observed, 
p(r) describes the dark matter density profile, and the 
integral is performed over the line-of-sight. dN 1 /dE 1 is 
the spectrum of prompt gamma rays generated per an- 
nihilation, which depends on the dominant annihilation 
channel(s). Note that Eq. [T] provides us with predictions 
for both the distribution of photons as a function of en- 
ergy, and as a function of the angle observed. It is this 
powerful combination of signatures that we will use to 
identify and separate dark matter annihilation products 
from astrophysical backgrounds [1] . 

With a perfect gamma ray detector, the distribution of 
dark matter events observed would follow precisely that 
described in Eq. Q] The LAT of the FGST has a finite 
point spread function, however, which will distort the 
observed angular distribution. In our analysis, we have 
modeled the point spread function of the FGST's LAT 
according to the performance described in Ref. Q. 

In addition to any gamma rays from dark matter anni- 
hilations coming from the region of the Galactic Center, 
significant astrophysical backgrounds are known to exist. 
In particular, HESS [3 and other ground-based gamma 
ray telescopes (TlJ have detected a rather bright gamma 
ray source coincident with the dynamical center of our 
galaxy (I = -0.055°, b = 0.0442°). The spectrum of 
this source has been measured to be a power-law of the 
form dNj/dEy « 10~ 8 GeV -1 cm -2 s _1 (£:/GeV)~ 2 ' 25 
between approximately 160 GeV and 20 TeV. Although 
HESS and other ground-based telescopes cannot easily 
measure the spectrum of this source at lower energies, 
it is likely that it will extend well into the range stud- 
ied by the FGST [12]], where it will provide a significant 
background for dark matter searches [l3j . Furthermore, 
gamma ray emission from numerous faint point sources 
and/or truly diffuse sources can provide a formidable 
background in the region of the Galactic Center, espe- 
cially near the disk of the Milky Way. To model this back- 
ground, we have studied the angular distribution of pho- 
tons observed by the FGST in the region of 3° < |Z| < 6°, 
and found that the emission is fairly well described by a 
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function which falls off exponentially away from the disk 
with a scale of roughly 2.2° to 1.2° for photons between 
300 MeV and 30 GeV, respectively. In addition to the 
previously described HESS point source, we will use this 
angular distribution as a template to model the diffuse 
background in our analysis. 

To start, we consider dark matter distributed accord- 
ing to a Navarro-Frenk- White (NFW) [1J] halo profile 
with a scale radius of 20 kpc and normalized such that the 
dark matter density at the location of the Solar System is 
equal to its value inferred by observations [15j . This pro- 
file, however, combined with the background model de- 
scribed above, leads to a distribution of photons that falls 
off more slowly as a function of angle from the Galaxy's 
dynamical center than is observed by the FGST. If we 
steepen the halo profile slightly, such that p(r) oc 1/r 7 , 
with 7 = 1.1 (7 = 1.0 for an NFW profile), the angu- 
lar distribution of the events matches the observations of 
the FGST very well. Such a steepening of the inner cusp 
could result from, for example, the adiabatic contraction 
of the profile due to the dynamics of the baryonic gas 
which dominate the gravitational potential of the Inner 
Galaxy 0. 

In Fig. [TJ we show the angular distribution of gamma 
rays as a function of angle from the Galaxy's dynami- 
cal center, as measured by the FGST over the period of 
August 4, 2008 to October 5, 2009, corresponding to an 
exposure in the Galactic Center region of approximately 
2.2 x 10 10 cm 2 sec at 300 MeV, and increasing to be- 
tween 3.5-4.3 xlO 10 cm 2 sec between 1 and 300 GeV. 
We include only those events in the "diffuse" class (as 
defined by the FGST collaboration) and do not include 
events with zenith angle greater than 105°, or from the 
region of the South Atlantic Anomaly. Below 10 GeV, 
we show separately those events that were converted in 
the front (thin) and back (thick) regions of the detector 
(these are treated separately to account for the differ- 
ing point spread functions for these event classes). In 
each frame, we compare these measurements to the an- 
gular distribution predicted for photon from annihilating 
dark matter with 7 = 1.1, the HESS point source, the 
diffuse background, and for the sum of these contribu- 
tions. In each case, these predictions take into account 
the point spread function of the FGST. For each range 
of energies shown, we have normalized the diffuse back- 
ground and dark matter contributions to provide the best 
fit to the data. Above 1 GeV, the normalization of the 
HESS source was determined by extrapolating the mea- 
sured power law spectrum to lower energies, which ap- 
pears to be consistent with the distribution observed by 
the FGST. Below 1 GeV, the FGST data do not appear 
to contain significant emission from this point source, so 
we suppress this contribution in these energy bins un- 
der the assumption that the extrapolation of the HESS 
spectrum becomes invalid around ~ 1 GeV. 

Thus far, we have performed these fits while remain- 
ing agnostic about the mass, annihilation cross section, 
and dominant annihilation channels of the dark matter 



particle, simply normalizing the angular shape predicted 
by our selected halo profile to the data. By comparing 
the relative normalizations required in each energy range 
shown in Fig. [TJ however, one can infer the shape of the 
gamma ray spectrum from our dark matter component 
and for the diffuse background. In particular, the rela- 
tive normalizations required of the diffuse component in 
the various energy ranges shown in Fig. [TJ imply that this 
background has a spectrum very approximately of the 
form dN 1 /dE 1 oc E~ 2 - 3 . 

In Fig. [2j we show the spectrum of photons observed 
by the FGST within 0.5° and 3° of the Galaxy's dy- 
namical center, and compare this to the spectrum in our 
best fit annihilating dark matter (plus diffuse and HESS 
point source backgrounds) model, including the effect of 
the FGST's point spread function (which noticeably sup- 
presses the lowest energy emission in the left frame) . The 
distinctive bump-like feature observed at ~l-5 GeV is 
easily accommodate by a fairly light dark matter parti- 
cle (rndm « 25 — 30 GeV) which annihilates to a bb final 
state. The normalization was best fit to an annihilation 
cross section of av w 9 x 10~ 26 cm 3 /s, which is a factor of 
about three times larger than that predicted for a simple 
s-wave annihilating thermal relic. 

It is interesting to note that the annihilation rate and 
halo profile shape found to best accommodate the FGST 
data here is very similar to that required to produce 
the excess synchrotron emission known as the "WMAP 
haze" 17]. To be produced by a dark matter particle 
as light as that described in this scenario, however, the 
observed hardness of the haze spectrum requires a fairly 
strong magnetic field in the region of the Galactic Cen- 
ter. It had been previously recognized that if the WMAP 
haze is the product of dark matter annihilations, then 
the FGST would likely be capable of identifying the cor- 
responding gamma ray signal [l8l ]. 

The low mass and relatively large annihilation cross 
section required in this scenario are also encouraging for 
the prospects of other gamma ray searches for dark mat- 
ter annihilation products, such as those observing dwarf 
galaxies and efforts to detect nearby subhalos. 

In conclusion, we have studied the angular distribu- 
tion and energy spectrum of gamma rays measured by 
the Fermi Gamma Ray Space Telescope in the region sur- 
rounding the Galactic Center, and find that this data is 
well described by a scenario in which a 25-30 GeV dark 
matter particle, distributed with a halo profile slightly 
steeper than NFW (7 = 1.1), is annihilating with a cross 
section within a factor of a few of the value predicted for 
a thermal relic. 

It should be noted, however, that if astrophysical back- 
grounds exist with a similar spectral shape and morphol- 
ogy to those predicted for annihilating dark matter (a 
spectrum peaking at ~ 1 — 3 GeV, distributed with ap- 
proximate spherical symmetry around the Galactic Cen- 
ter proportional to r~ 2 2 ), the analysis performed here 
would not differentiate the resulting background from 
dark matter annihilation products. Gamma rays from 
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FIG. 1: The angular distribution of gamma rays around the Galactic Center observed by the FGST. In each frame, the dashed 
line denotes the shape predicted for the annihilation products of dark matter distributed according to a halo profile which is 
slightly cuspier than NFW (7 = 1.1). The dotted line is the prediction from the previously discovered TeV point source located 
at the Milky Way's dynamical center, while the dot-dashed line denotes the diffuse background described in the text (which, 
although included in each case, falls below the range of rates shown in the upper four frames). The solid line is the sum of 
these contributions. 
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FIG. 2: The gamma ray spectrum measured by the FGST within 0.5° (left) and 3° (right) of the Milky Way's dynamical 
center. In each frame, the dashed line denotes the predicted spectrum from a 28 GeV dark matter particle annihilating to 
bb with a cross section of av = 9 x 10 -26 cm 3 /s, and distributed according to a halo profile slightly more cusped than NFW 
(7 = 1.1). The dotted and dot-dashed lines denote the contributions from the previously discovered TeV point source located 
at the Milky Way's dynamical center and the diffuse background, respectively. The solid line is the sum of these contributions. 



pion decay taking place with a roughly spherically sym- 
metric distribution around the Galactic Center, for ex- 
ample, could be difficult to distinguish. Further informa- 
tion will thus be required to determine the origin of these 
photons. 
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